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Abstract 

We report new results for the Thomas- Fermi groundstate of a spin- polarized dipolar interacting 
Bose-Einstein condensate for the case when the external magnetic field B is not orientated parallel 
to a principal axis, but is aligned parallel to a symmetry plane of a harmonic anisotropic trap. For 
a dipole interaction strength parameter eo 7^ the release energy of the condensate depends on 
the trap orientation angle "dx between the principal axis e^^T of the trap and the field B. From 
the quasiclassical Josephson equation of macroscopic quantum physics we determine the low-lying 
eigenfrequencies of small amplitude collective modes of the condensate density for various trap 
frequencies Ua and trap orientation angles 'dx- For the special case of a spherical harmonic trap with 
trap frequency uj it is rigorously shown for — i < e/) < 1, that a pure s-wave symmetry breather 
excitation of the condensate density exists, that oscillates at a constant frequency Vis = V^oj 
around the groundstate cloud, despite the well known fact, that the shape of the groundstate 
cloud of a spin-polarized dipolar condensate is for ed 7^ not isotropic. For i?t 7^ the small 
amplitude modes of the particle density with isotropic and quadrupolar symmetry consist of two 
groups. There exist four modes that are combinations of basis functions with s-wave, dj.2_y2 - and 
-wave, and also daj^-wave symmetry, and two modes that are combinations of basis functions 
with dyz- and dxy-wave symmetry. A characteristic diff^erence in the dependence of the frequencies 
of these six collective modes on the dipole interaction strength parameter ed for prolate and 
oblate harmonic tri-axial traps, respectively, is suggested to be used as an experimental method to 
measure the s-wave scattering length of the atoms. 
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I. INTRODUCTION 



Experiments with trapped, extremely dilute gas clouds, consisting of identical atoms 
with mass m*, and forming at ultracold temperatures a quantum degenerate Bose-Einstein 
condensate (BEC), are nowadays a research focus in many laboratories. Early on it has 
been realized that cold atom clouds do not form an ideal Bose gas, but experience in the 
low energy sector of the system isotropic interaction forces that can be well described by 
a microscopic s-wave scattering length While the size of an ideal Bose gas confined 



inside a harmonic trap with trap frequency u is determined by the width = y 
of the groundstate wavefunction of a single particle, the size of an interacting cold atom 
cloud consisting of a large number >> 1 of condensed Bose atoms may increase to much 
larger distances A^-^ = l^]rNas_y_ Fortunately, the necessary requirement ^ 1 

for observing a BEC in a harmonic trap can be realized simultaneously with the condition 
of a small diluteness parameter tt-qO^ <^ 1, so that the mean field theory of Ginzburg and 
Pitaevskii for interacting Bose systems is applicable for a wide range of parameters a^^ and 
Os- As the length Ktf increases with increasing A^ the kinetic energy Ek — 2n^\^ 
interacting particles in the groundstate eventually becomes much smaller than the potential 
energy Vt — ^u^A'^p of the particles, because the density inside the BEC becomes a 
smooth and slowly varying function of position. In the Thomas-Fermi approximation the 
kinetic energy term for the particles in the groundstate of the BEC is neglected altogether. 
This is justified when the chemical potential /i of the interacting system is much larger 
than the chemical potential ~ of the non interacting Bose gas. So for ^^^"^ ^ 1 the 
dominant balance required for mechanical equilibrium of a trapped BEC is between the 
repulsive interactions of the atoms and the confinement forces of the trap. 

New interesting physics can be observed when in addition to the usual s-wave contact 
interaction the atoms get influenced by long ranged dipole-dipole forces Sj. This occurs, for 
example, for Bose atoms with nuclear spin / = and integer (electronic) spin S, thus giving 
rise to a multiplet —S < Ms < S of atomic magnetic dipole moments with ^-component 
2^bMs. a transition metal atom like chromium ^^Cr has / = and S = 3. On the other 
hand. Alkali atoms like ^'^Rb carry (nuclear) spin / = | and 5" = |, thus coupling to a total 
spin F = Hn the lowest energy state. The first experimental stndy ot ntagnetrc dipole-dipole 
interactions in a BEC was realized with Cr- atoms |4|] carrying a large magnetic moment 



I (M) I = 6/iR. A quantum degenerate F = 1 spinor BEC was synthezised successfully with 
^^i?6-atoms(5|. Recently, intrinsically anisotropic BEC systems with electric dipole-dipole 



interactions between polar molecules have been studied experimentally [20[. New research 
directions are concerned with magnetic quantum gases consisting of heavy rare earth atoms 
like Thulium^] with |(M)| = 4//^ , Erbium|9;] with |(M)| = T/x^, and Dysprosium [Sj] with 
|(M)| = lO/ifi. 

In the ensuing considerations we study spin-polarized Bose atom clouds. When the 
magnetic dipole moments of the atoms are 100% polarized under a homogeneous external 
magnetic induction field B = B^'^^^^e^ , so that all atoms in the cold gas cloud carry the 
identical effective magnetic moment, say (M) = —2^bS , it is still possible to describe 
the Bose condensed groundstate ^ oi N interacting atoms by a scalar Hartree ansatz 

^ (r«, .., r(^)) = ^(r(i)) ■ ^(r^^)) . . . ^(r^^)) (1) 

The expectation value of the many body Hamiltonian H , evaluated with such a trial wave 
function \1/ consisting of a product of identical one-particle wave functions ip{r), is then 
minimized with respect to variations of that one-particle wave function iI){y). The optimal 
one-particle wavefunction ip{v) so found is a solution to the Gross-Pitaevskii equation 



(2) 



-Vl + Vrir)- ^i+{N -I) ! dV[/(r,r') 



/m2 



|^(r')| 



^(r) = 



2m* 

Here, VT(r) denotes the potential of the trap, and U (r, r') describes the interaction potential 
between two Bosons. The chemical potential /i is a Lagrange parameter connected to the 
particle number in the condensate by the constraint: 

d?r' \ip{r'f = 1 (3) 

II. THOMAS-FERMI THEORY OF SPIN-POLARIZED DIPOLAR BOSE-EINSTEIN 
CONDENSATE 

In the following we investigate the macroscopic quantum degenerate groundstate of a 
spin-polarized system of interacting Bose atoms carrying a magnetic dipole moment |(M)|. 
The interaction potential 

[/(r,r') = f/o(r,r') + [/^d(r,r') (4) 



between two atoms, one at position r and the other at r' , consists of two contributions, the 
short ranged isotropic s— wave interaction pseudopotential 



[/o(r,r')= (7. (r - r') 



(5) 



Qs = a« 



, and the long ranged magnetic dipole-dipole interaction potential: 



Umd (r, r ) 



9md 



4tt 

9md = 1^0 I (M) I 



r — r 
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./|3 



3 {r, - 

I /|5 

r — r 



(6) 



Here the external magnetic induction field B is orientated parallel to the Cartesian unit 
vector in the laboratory frame, so that the magnetic moments of two interacting atoms, 
one at positions r and the other at r', are both aligned parallel to e^. 
Using well known identities 



r — r 



/|3 Ij. 



1 



(7) 



-V 



r — r 



r — r 



and introducing the dimensionless parameter 



Qmd _ ^0 |(M)|' 



(8) 



as a measure of relative strength of magnetic dipole interaction forces, the interaction po 
tential between two atoms in the gas cloud may be rewritten in the guise: 

1 1 



9s 



3e 



D- 



(9) 



dr"^ An |r — r' 

In the Thomas- Fermi approximation the particle density profile in the groundstate of the 
trapped BEC 



n r 



Nip [r] 



is a solution to the integral equation 

1 



(1 - Ed) riTF (r) - Sed 



drl An 



r — V 



-riTF (r 



(10) 



(11) 



This is actually a non linear problem, because the solution of the integral equation is sought 
inside the Thomas-Fermi cloud 

Btf = {r eR^\ riTFir) > O} (12) 

, which region is not known a priori. The determination of the shape of the cloud , or 
its boundary c?Dtf , is part of the problem. 

Eberlein et al. have found for a dipolar interacting BEC confined inside a harmonic 
trap, that despite the non local anisotropic dipole-dipole interaction term, the domain Htf 
always maintains the shape of an ellipsoid, as in the case Ed = 0, but with different semi axes. 
We confirm this finding and present new results for the case, when the external magnetic 
field B is not in alignment with the principal axis e^^r of the trap. 

Consider a harmonic anisotropic trap potential Vrir) with its minimum at position r = 
, and with the principal axis e^^r of the trap not in alignment with the field B : 

Vrir) = ^ {ulr^T + ^VIt + r) (13) 

For Ux ^ ojy , ujy ^ ojz and Uz 7^ uJx surfaces of constant trap potential VT(r) = ^ > 
have the geometrical shape of a tri-axial ellipsoid. Three mutually orthogonal Cartesian 
unit vectors e^^x , Gy,T and e^^^ determine the orientation of the principal axes of such a 
trap. The magnetic field B is then in general a linear combination of all three principal axis 
vectors: B =Brc^T^x,T+ By^TGy,T + Bz^t^z,t- For simplicity we restrict our considerations in 
the following to the special case, when the magnetic field B and the principal axis ez,T of 
the trap span a symmetry plane of the trap, say the plane = 0. Then we have in ( !T3l) : 

(14) 

rx,T = Tx (^?t) = COS {dr) Tx + sin {^t) Tz 
ry,T = Ty {dr) = ry 

rz,T = (^?t) = - sin (dr) rx + cos (dr) rz 

, i.e. the principal axis e^^r of the trap is turned by an angle {}t around the rotation axes 
ey^T -L B (see Fig.l). 

As is indicated in Fig.l, the selfconsistent solution of ( ITT]) for the density distribution 
utf (r) reveals, that the principal axis e^^o of the Thomas- Fermi cloud Dtf is rotated away 
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FIG. 1. (Color online) Orientation of principal axis e^.T of harmonic trap and orientation of 
principal axis e^^o of Thomas-Fermi ellipsoid Dtp relative to the spin polarizing external magnetic 
field B. The inset on the left hand side corresponds to a cut of Otf with the symmetry plane 
y = 0. Arrows indicate the orientation of the spin polarizing magnetic field B relative to the 
principal axis e^^o of Dj-f- 

from the direction of the external field by an angle i^o 7^ ^t- Accordingly, the density profile 
associated with the ellipsoidal domain has the general form 

-^x \ J 

where 

(16) 

rx = r^ = cos {'do) + sin (^?o) rz 
Ty {^o) = ry = ry 

rz (i^o) = rz = - sin (^9o) + cos (^9o) ^z 

Only in the highly symmetric case "i^T = the principal axis vector e^^o of the ellipsoid 3tf 
is orientated parallel to B. For < t?T < f it is found from the selfconsistent solution for 
the density profile 77,^^(1"), that ^ ■, i.e. the principal-axis e^^o of '^tf is never in 
alignment with the field B, nor is it in alignment with the principal axis ^ of the trap 
(see Fig.l and Fig. 2). 
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FIG. 2. (Color online) Angle difference 'do — t^t vs. trap orientation angle 'St for selfconsistent 
groundstate density profile utf (r). The curves shown correspond to a dipole interaction strength 
£d = 0.5 and to various frequency ratios u^x '■ '■ of the harmonic trap: a) 2 : 6 : 3 purple 
empty diamond ; aa) 3:6:2 purple full diamond; b) 3 : 2 : 6 red empty square; bb) 6:2:3 
red full square; c) 1 : 2 : 2 blue empty circle; cc) 2:2:1 blue full circle; d) 2 : 3 : 6 green empty 
triangle; dd) 6 : 3 : 2 green full triangle. 



The particle density distribution nrir(r) inside the Thomas- Fermi ellipsoid is strat- 
ified. Like an onion it consists of a series of thin homoeoidal shells of constant density 



nTpij) = '"'0 (l ~ i^^) = const 
< < 1 



(17) 



Strata of equal density thus correspond to ellipsoidal shells concentric and similar to the 
bounding ellipsoidal shell 9Dtf , but with scaled semi- axes z/Aq. 
For a general tri- axial ellipsoid J^tf the normalization integral 



(fr'riTF (r') = 



leads to 
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(18) 



(19) 



So, the problem is to determine from ( ITTj) the three semi axes A^; , Ay , , the orientational 
angle 'do and the chemical potential /i. 

As has been emphasized by Eberlein et al. in ^ the central task in solving the Thomas- 
Fermi integral equation flTT]) is to calculate the potential function 

0TF (r) = ^ / d^r'^^TiTF (r') (20) 
4vr Jb^^ \r - r'\ 

for a heterogeneous particle density distribution utf (r). Then, because 

- V^(Ptf (r) = riTF (r) (21) 
, the Thomas-Fermi integral equation ( ITTl) 



(1 - Sd) utf (r) - 3eD^<PTF (r) = ^^-^ (22) 
becomes (in free space) equivalent to a partial differential equation of potential theory: 

(23) 



9s 

In order that the differential operator on the left hand side is positive definite it is required 
that 

-l<en<l (24) 

Due to the (i22-anisotropy of the dipole-dipole interaction between two spin-polarized atoms 
at position r and r' the dipole-dipole interaction part of the potential U (r, r') is attractive 
or repulsive, depending on the orientation of the distance vector r — r' relative to the field 
vector B. So, if ££> > 1 or ed < — | , attractive forces prevail and the system will collapse. 
Of course, a better criterion for stability is to calculate the frequencies of the collective 
modes of the dipolar interacting BEG. We shall present in the next section new results for 
quadrupolar like modes of the density fluctuations around the groundstate density profile 
Utf (r)- 

The differential equation (123!) makes it manifest, that the integral operator with in- 
tegration domain B)tf and kernel in the Thomas- Fermi integral equation fl22] ) for 
r gDtf maps a quadratic form htf (r) spanned by the linearly independent basis functions 

, T^ai^b} a<b<e{x y z} ^^^o a quartic form (pTF (r) spanned by linearly independent basis func- 
tions {1, rarj,, rarbrcrd}a<b<c<de{x y z}- Because for a harmonic trap potential Vr (r) the right 
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hand side of (123|) is (by definition) a quadratic form, the problem would be exactly solved, 
provided the coefficients of the quadratic form presented by the derivatives J^-^tf (r) of the 
potential function ( l20l) can be found. 

We describe now a very convenient method to determine the coefficients of the quadratic 
form ■^4'TF (r), which is all we need to solve the Thomas- Fermi integral equation (1221 ). As 
a matter of fact, three-dimensional integrals of the type 



47r In s - 



s' 



n'2 



-,12 



A? 



(25) 



/ = 0,1,2,3,... 



s' G 



J2 o'2 „/2 
a; y z 



have been calculated analytically by S. Chandrasekhar 15| in his magisterial treatment of 
the ellipsoidal figures of equilibrium of gravitating and rotating gas clouds in astrophysics. 
He showed, that (s) can be represented exactly in terms of singularity free fast convergent 
one- dimensional integrals. Chandrasekhar's result for the three-dimensional integral $/ (s) 
at an internal point of the ellipsoid D is: 



s e Btf 

I = 0,1,2,3,... 



(26) 



(s) 



Xx^y^z 1 

4 l + l 



du 



' ^Jm + u) {XI + u) {Xl + u)] 



XI + u Xl + u Xl + u 



i+i 



To solve the Thomas-Fermi integral equation we now make use of this result for the 
special case 1 = 1. Because the Cartesian coordinates of a point presented in the principal 
axes frame of the Thomas- Fermi ellipsoid Dtf are connected to the Cartesian coordinates 
ra of that same point in the laboratory frame by a rotation. 



Tb (^?o) = n 



Gy) Ta 

a(i{x,y,z} 



(27) 



Tlba (^?o; % 



COS ("i^o) sin (-i^o) 

1 
— sin (-i^o) cos {do) 



J ba 
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, and taking into account that under such a rotation TZ {-^o', ^y) we have |r — r'| = |r — r'|, 
we immediately see that the potential function 

(/)TF(r) = no$i[7(r)] (28) 

at a position r G Dtf is a quartic form with regard to the linearly independent basis 
functions {1, '^a^'fe, ^a'"fe'"c?^d}a<b<c<de{a; y z}- ^ sccoud ordcr derivative of a quartic the 
function ^0tf (r) is then manifestly a quadratic form, spanned by a linear combination of 
the basis functions {l, r^, r^, rf, rj-r^} , or taking into account (1271) . it is spanned by a linear 
combination of the basis functions |l, r^, r^, r^, r^jr^}. The transformation from one basis 
system to the other is accomplished by the orthogonal transformation (!27|) . 

The coefficients Cab of the quadratic form 

reDTF (29) 

d . . TIq / 2 2 2 \ 

-q:^9TF (r) = — (-Coo + C^xT^ + CyyTy + C^^r^ + C^^zTxTz) 

depend on the trap orientation angle i^o and the semi- axes A^: , , A^ of the ellipsoid D^f 
via the following one- dimensional integrals: 

a, 6 G {x, y, z] (30) 



-^a (Aj^, Aj;, A^) — \x^y^z \ 

Jo 



du 



-^afe (Ax, Ay, A^) = A^Aj/A^ / 

Jo 



du 



^i>^l + u) {XI + u) {XI + u) + + 



In the appendix |X] some of the properties of these so called index integrals are listed. We 
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find 



(31) 



Coo = sin^ (^9o) 4 + cos^ {^q) 

cos2 i^o) sin^ (4) (4,. + h,) + [cos^ (4) + sin^ (^9o)] /.x 
+2 sin^ (^9o) cos2 (^9o) (4. - 2/,, + /,,) 

Cyy = sin^ (t9o) 4y + cos^ {^o) hy 

c^^ = 3 ■ [sin^ {'do) 4, + cos^ {do) l^^ + 2 sin^ {^o) cos^ (t^o) 4x] 
Ca;z = 6 sin {do) cos {do) [sin^ (^9o) (4x - ^zx) + cos^ (^^o) (-^^x - lzz)\ 



It is advantageous to work in the geometry under consideration not with the basis func- 



tion of the coordinate system around the axis e.y by the trap orientation angle "i^o as defined 
in ( !27l) . The exact solution of the Thomas- Fermi integral equation ( l22|) is then obtained in- 
serting the corresponding explicit expressions for the quadratic form ^0tf (r) and the trap 
potential VT(r). From the condition, that the prefactors of the linearly independent basis 



tions {I,r2,r2,r2,r,r^ 





r^r^} in fl22|) should vanish identically, the following set of coupled 




(32) 



(33) 




2^ 



{ul cos^ {dT - do) + sin2 {dT - do)\ 
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I 



(34) 



+ 3|o [cos^ (^?o) /.J, + sin^ (??o) 4s,] 



(35) 



> no = — [u;^ sin^ (^?t - A) 
^9s 



+ 3|D [3cos2 (^o)/z. + sin2 ^^^^^^^j 



UjI COS^ (l?T - l?o)] 



(36) 



3£d 



* 2 2 

Sin (2^o) hz no = ^ "^"""^^ • sin (2^t - 2^?o) 



2gs 



The normalization condition connects the density Uq at the center of the Thomas- Fermi 
domain D^^p to the product of the semi-axes: 



no = 



15 



(37) 



Let us first write the selfconsistency equations without dipole interaction setting Ed — 0. 
There follows 



n, 



(0) _ 



(0) 



m 



Using the normalization 



and introducing the definitions 



(0) 



CO. 



nn = 



2//(o) " 



15 



(38) 



(39) 



h 



rrru 



(40) 
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we obtain for the chemical potential /i'^^-' and the semi-axes X^a ^ of a BEC inside an anisotropic 
harmonic trap in the Thomas-Fermi regime well known results: 



Ai°UA^)=(^^) =^A (42) 



A = ttn 



15 47r iVa,\ ^ 
47r a,„ 



The exponent | is characteristic for the large N scaling of the chemical potential n^'^^ (N) 
of a BEC confined inside a harmonic trap |l|. 

The ensuing calculations simplify making use of elementary scaling relations that hold 
for single index and double index integrals: 

la (A., A„ A,) = 4 ^, 1^ = 7, (43) 

Ac lab (Ax, Aj,, Xz) = lab ''"^ ^ )^ '^"^ 



Using the obvious relation 



uq aI ''Ay ■'aI 



nl°^ Aa^-Aj^A^ 



(44) 



the selfconsistency problem posed by fl32l) - fl36l) may then be reduced to three coupled equa- 
tions for the ratios , and the equilibrium orientation angle -i^o as functions of the trap 
orientation angle , the trap frequencies uja and the dipole interaction strength parameter 
Sd : 

(45) 

XI _ ul sin^ {dT - 4) + ujI cos^ {dT - ^o) ^-^D + [cos^ (^?o) + 3 sin^ (^p) 

XI ~ ul cos2 - ^o) + sin^ - ^9o) ' 1 - £d + ^ [3 cos^ (^o) + sin^ {d^) 7,,] 



(46) 

XI _ ul sin^ (^T - ^o) + cos^ (t^T - ^o) 1 - + ^§ [cos^ (^o) 7,, + sin^ (t^q) 7,,] 
A^ ~ a;2 ■ 1 _ + ^^^2 (^^) j^^ + gj^s (^^) J^^] 



13 



(47) 



tan (2i?o 



[ujz 



sin (2t9' 



T 



l-eD + 



■1 \^ 



COs2(.^lo)/.y+sin2(i?o)/^y 



We have found, that the set of selfconsistency equations fHSj) . fH6|) and P7|) may be conve- 
niently solved numerically by the method of fixed point iteration. Using identities like 

1 1 1 

cos^ (i?o 



2 2 ^1 + tan2 (2i9o) 

the evaluation of trigonometric functions in the iteration process can be completely avoided. 

Once the ratios y' i ^ the orientation angle ?9o are known, it follows directly from 
(1351) and (EHl): 



(48) 



sd + 



Set 



[3cos2 (4)/,, + sin2 (^q)/^,^] 



l?o) + C0S2 - 1?o) 



A2 A5; 



For the special case of an isotropic harmonic trap [l| the principal effect of the dipole- 
dipole interaction on the groundstate density profile utf (r) of a dipolar interacting BEG 
is the well known elongation of the semi- axis parallel to B, and the distortion of the 
semi-axes A^^ = A^ perpendicular to B towards smaller values: 



Ax _ 1 — l^D + ••■ 

aT ~ i + hD + ... 



(49) 



5^ 

In Fig. 2 and Fig. 3 self consistent solutions of the coupled equations (H5|) . ( H6|) and ( H7|) for the 
equihbrium angle ^q—'^t and the semi-axes Aa are plotted as functions of the trap orientation 
angle for various magnetic dipole interaction strength parameters Sd assuming a tri-axial 
trap anisotropy ratio Ux '■ ojy : Uz = Q : 3 : 2 . 

Once the semi-axes Aa of the particle density htf (r) ( pTSl) are determined, then (jS]) 
gives us the value uq of the density of the ellipsoidal shaped BEG at its center. In Fig.4 
the ratio is plotted vs. the trap orientation angle i^t , the inset showing for selected 
trap orientation angles t?T ^ {O' f' f } ^^^^ of the corresponding selfconsistently determined 
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FIG. 3. (Color online) Plot of semi-axes A^; , Aj^ and A^ of selfconsistent groundstate density profile 
n^Flr) vs. trap orientation angle for various dipolar interaction strength ££>: a) en = 0.2 
dotted line ; h) £d = 0.5 dashed line; c) £d = 0.8 solid line. The ratio of trap frequencies is 
Ux : ijJy '■ U!z = G ■ 'i '■ ■ 

Thomas- Fermi ellipsoid Dj-f with the symmetry plane y = 0. For a prolate trap the density 
is largest for i!}t = , because then the net mutual dipole force between atom pairs inside 
the domain H^f is attractive. As i!}t increases the density no becomes smaller and assumes a 
minimum at "i^y = 90° , because for a parallel alignment the net mutual dipole force between 
atom pairs inside the domain 3tf is repulsive. 

Finally, there follows from (!32|) an explicit formula for the chemical potential: 

fi=ll-en + lsD [sin' (^o) h + cos' {^o) 7.] | ^ ■ f^^'^ (N) (50) 

The dependence of fi on particle number is solely described by the factor (N) , i.e. 
the ratio jjgj is independent on particle number A^. In Fig. 5 the chemical potential fj, 
is plotted vs. the trap orientation angle for different values of the dipole interaction 
strength e^^. While for an isotropic harmonic trap the chemical potential /i doesn't change 
to first order in Sd , one finds for an anisotropic harmonic trap, making a straightforward 
expansion to the first order in the dipole interaction strength Ed , for the case of an oblate 
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1.25 



sd = 0.5 




FIG. 4. (Color online) Selfconsistent particle density uq at center position of anisotropic harmonic 
trap plotted vs. trap orientation angle i?^. The inset displays cuts of the Thomas-Fermi ellipsoid 
^])TF with the symmetry plane y = for trap orientation angles 'dx = 0°, 'dx = 45° and = 90°. 
The ratio of trap frequencies is ojx ojy : ojz = 2 : 2 : 1 ^ the dipole interaction strength \s ed = 0.5 



(pancake shaped) trap that 



deo 



£15=0 



> 0, and for a prolate (cigar shaped) trap that 



d^ 

deo 



£13=0 



< , respectively. 



The displayed characteristic dependence of chemical potential /i on the trap orientation 
angle dj- should be observable as the release energy of a spin-polarized dipolar interacting 
BEC confined in a harmonic trap, when the trap potential is suddenly switched off to zero, 
and subsequently the dilute atom gas cloud undergoes a ballistic expansion We find 
within the range of validity of the Thomas-Fermi approximation that there holds also in the 
presence of long ranged dipole-dipole interactions 



Er = E; 



int 



(51) 
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FIG. 5. (Color online) Chemical potential /x vs. trap orientation angle ■i?T' for different values of 
dipole interaction strength, a) ed = 0.2 red square ; b) = 0.5 blue circle ; c) = 0.8 purple 
triangle. The ratio of trap frequencies is cj^ : : = 6 : 3 : 2 . 

As a matter of fact, the total energy 

E = (^^^^ = (^Hkin + Hpot + Hint^^ (52) 

, and the interaction energy Eint = (^HintJ are connected in the groundstate ^ of the 
BEC, see ([H), by the general relation 



E = fiN- Eint (53) 

This applies, because the optimal one-particle wave function if) (r) building the particle 
groundstate ^ solves the GP equation, so that the expectation value (^Hkin^ of the kinetic 
energy can be re-expressed via the GP equation in terms of the chemical potential /i and 
the interaction energy Eint ■ On the other hand, the total energy E of the BEC is connected 
to the chemical potential by the general relation 

It follows from (150!) . and the established scaling ( BTj) of (A^) oc Ns for a large particle 
number A^ ^ 1, that up to a constant that is independent on A^ there also holds for a 
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spin-polarized dipolar BEC confined inside a harmonic trap the well known relation 



E = ^/iiV (55) 

, and therefore: 

E,^t = fiN-E = ^fiN (56) 



III. COLLECTIVE MODES OF SMALL AMPLITUDE DENSITY OSCILLATIONS 

A. Parametrization of Low-Lying Excitations in Tri- Axial Harmonic Trap 

An important test of the macroscopic quantum physics of a BEC is the study of ele- 
mentary excitations above the groundstate. One technique to excite low energy collective 
modes of a BEC is to suddenly modify the trap potential. For example, shifting the center 
of the trap excites the dipole modes, that is the motion of the center of mass of a BEC cloud 
around its equilibrium position in a harmonic trap. Changing the curvature of the trap by 
switching the trap frequencies may excite the breather mode. In an anisotropic harmonic 

n 

trap there also exist the so called sissors modes [12|, which can be excited by rotating a 
principal axis of the trap, thus pushing the atom cloud in the trap away from equilibrium. 
A recently reported elegant new experimental technique excites a BEC by modulating the 
ield dependence of the atomic scattering length near to a magnetic Feshbach resonance 



3, 



lid 



We calculate in this section for the case of a harmonic anisotropic trap with arbitrary 
trap orientation angle 'dx the small amplitude collective modes of a dipolar interacting 
spin-polarized BEC at very low energy, so that the wavelength of the excitations becomes 
comparable to the size of the Thomas- Fermi length A. For a large number >> 1 of 
particles in the BEC the particle density n (r, t) and the macroscopic Josephson phase S (r, t) 
are conjugate variables, so that the collective dynamics of the system (ignoring a small 
quantum pressure) is governed by the standard canonical equations of motion of macroscopic 
quantum physics [ij: 

(57) 
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(58) 



— nr,t +— > — 



n{r,t) 



dS (r,t) 







The gradient of the phase, the velocity field v(r,t) = ^V(5s(r,t), is directly connected 
to the density of particle current, j (r, t) = n (r, t) v (r, t) [1]. The quasiclassical Josephson 
equation and the continuity equation need to be solved subject to the normalization 
condition 

/ d^r n (r, t) = N (59) 
, where for a fiuctuating particle density n (r, t) the positivity domain 



D(t) = {r G M^l n(r,t) > O} 



(60) 



depends ( in principle) on time too! 

Let us assume a perturbation expansion of phase and density of the form 



S{v,t) = -^-t + 5s (r,t) 

h 



(61) 



n (r, t) = riTF (r) + (r, t) 

, with 6s (r, t) and 5n (r, t) denoting small fluctuations of phase and particle density around 
the groundstate of the BEC. 

For convenience we work from now on, unless otherwise explicitely stated, in the prin- 
cipal axes frame of the Thomas- Fermi ellipsoid V>tf- As is indicated in Fig. 1, for a trap 
orientation angle ?9t 7^ , the magnetic field B is not aligned parallel to the principal 
axis vector e^^o of the ellipsoid 3tf- In this case the interaction energy U{r,r') between 
two spin-polarized atoms at position r and r' , both carrying a magnetic dipole moment 
(M) = {2hbS) m orientated (anti-)parallel to B is then given by 



U (r, r') = 



ez?) <^('Ur - r') - Se^ (m ■ V^^' ^ ^ 



(62) 



47r |r — r' 

In the geometry under consideration the axis q of the ellipsoid "^tf is rotated around the 
axis Gy^T of the trap by an angle do , and we assume m ± : 

f)2 f)2 C)2 

(m ■ Vr)' = sin^ (4) — + cos^ (^^o) ^ + sin (2^o) (63) 

19 



Upon linearization of ( !57|) and (158|) we reproduce to order zero in the expansion the Thomas- 
Fermi integral equation f|TT]) determining the equilibrium density profile utf (r) and the 
chemical potential discussed already in the previous section. 

To derive the equations of motion for the fluctuations of the phase 5s (r, t) and the density 
5n (r, t) let us first consider for a kernel K (r, r') that couples to particle density n (r', t) the 
associated fluctuation 



5K{t)= f (fr' K{r,r')n{r',t) - I d^r' K {r,r') utf {r') 

Jn(t) Jdtf 



(64) 



}{t) JOtf 

In principle there are two contributions to 6K (t) . One is generated by the time dependence 
of the density distribution, the other results from a change of the integration domain D (t). 
Introducing the Heaviside distribution 

1 + sign(x) 



e^, (x) 



(65) 



we rewrite f l64|) as 



dK (t) = J d^r'K (r, r') < 



Qh [utf (r') + Sn (r', t)] [utf (r') + Sn (r', t)] 



-Qh [utf (r')] riTF (r') 
An expansion to the first order in the small quantitiy 6n leads to 



(66) 



6K (t) = j d^'K (r, r') 6n (r', t) {6 [utf (r)] + utf (r) 6d [utf (r)]} + o {\6nf) (67) 



where 6d (x) denotes the Dirac delta-distribution: 



d 



(68) 



Here, the term proportional to utf (i"') Sd [utf (r')] corresponds to a surface integral over 
the boundary OI^tf of the Tomas- Fermi domain IDtf- However, because utf (r') = for 
r' G d^TF the value of this surface integral is zero. This means the fluctuation of the 
integration domain © (t) around the shape of the equilibrium cloud 3tf as caused by a 
density fluctuation 6n = n (r, t)—nTF (r) represents only a small correction to 5K {t) beyond 
first order accuracy: 



6K{t) 



dVir(r,r') 5n {r' ,t) + o {\6nf) 



(69) 
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It follows now directly from (157|) for the time derivative of the phase fluctuation 5s (r, t) 
substituting into fl69l) the special case K {r,r') = U (r,r'), that 

d 



a r 

h—Ss (r, t) = ~ d?r' U (r, r ) 5n (r , t) 



and further from (l58 







a e {x,y,z} 

— 5n(r,t) + — — 



(70) 



(71) 



This is the well known wave equation describing the collective density excitations of a BEC in 
the quantum hydrodynamic limit [l| for the case of a non local interaction potential U (r, r'). 

In the ensuing discussion we consider the collective modes of small amplitude oscilla- 
tions around the equilibrium density distribution htf (r) in the trap. The corresponding 
density fluctuations 6n (r, t) may be expanded with respect to a set of linearly independent 
multinomials {1, r^, r^rb, ranVc, ...}„ 



<b<c.. 



a, 6, c e {x, ?/, z} (72) 
5n (r, t) = [t) + ^ (t) + ^ ^n^^ (t) r^n + ^^^Sl, (0 ranr,... 

a a,b a,b,c 

We restrict here to the terms of zero order , first order and second order of the density 
fluctuation 6n (r, t) as described by amplitudes 6n^^^ (t) , 6n^a ^ (t) and 6n^^j^ (t) , respectively. 
Third and higher order terms proportional to Sn^^^j^ (t) etc. we shall consider elsewhere. 

A very convenient approach to the determination of the small amplitude oscillations 
6n (r,t) around the equilibrium density Utf (r) = nxp (r; A) is to parametrize the density 
fluctuations 6n (r, t) in terms of a displacement vectorfield i] (r, t) and in terms of a dilatation 
amplitude vector ( (t) : 



6n (r, t) = Utf [r + (r, t) ; A + C {t)] - utf (r; A) 



(73) 



Neglecting small higher order terms regarding the size of the amplitudes {rjl and \(\ , a 
general density fluctuation 6n (r, t) around the groundstate of the BEC cloud is then 



b E {x,y,z} 

driTF (r; A) 



(74) 



5niv,t) = J2 



Vb (r, t) + — Cb [t) 



+ o{\v\' + \Cf) 
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It follows directly from atom number conservation, and substituting into (169|1 the special 
case K (r, r') = 1 , that 

[ d^r'6n{r',t) = (75) 

Upon insertion of f l71|) into fl7^ , and using the theorem of GauB, we see that the displace- 
ment vectorfield r] (r, t) is necessarily a solenoidal vectorfield: 



div?] (r, t) = 



(76) 



Making a partial integration this property of rj (r, t) implies for the phase fluctuation ([70 



(77) 



d 



4^^dVr^,(r',t)^n^Wr';A) 
-C^ it) A /b,, U (r, r') (r'; A) 



A simplification results for standard two body interaction forces among two atoms, say at 
position r and r' , that obey Newton's law "actio = reactio": 



dU (r, r') dU (r, r') 



dr'. 



dn 



(78) 



Then a general density fluctuation 5n (r, t) around equilibrium, as parametrized by fl74|) in 
terms of a solenoidal displacement rj {r,t) and a dilatation amplitude ( (t), is connected to 
the time derivative of the fluctuation 6s (r,t) of the Josephson phase by 



b e {x,y,z} 



(79) 



d 



i:L^d^r'Uir,r')nTFir';X) r^,{v\t) 



+Cb{t)i-Ld^r'U{v,v')nTF (r';A) 



Indeed, this exact representation of the quasiclassical hydrodynamic phase fluctuations of a 
BEG for interactions U (r, r') that obey (!78|) represents a convenient starting point for our 
analytic calculation of the collective density oscillations of a spin-polarized dipolar interact- 
ing BEG in a trap. 

Like the density fluctuations in (I72|) we may also expand the Gartesian components 
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rja (r, t) of the displacement vectorfield: 



a,b,ce {x,y,z} 



(80) 



Va (r, t) = r/f ) (t) + J2 Vat (i) n +^2'^. 



(t) nrc.... 



b b,c 



The zero order term ija (t) describes a homogeneous displacement of the center of mass of 
a BEC cloud in a trap. Actually, any system of particles that interact via two-body forces 
obeying to (!78|) . has the property that the motion of the center of mass separates from the 
equations of motion of the other degrees of freedom of the system. As a result the frequency 
of the dipole modes of a trapped atom gas cloud is independent on any such interactions, 
because the center of mass of the cloud moves like a single particle of mass Nm* in the 
external trap potential Vrir). For a harmonic trap the frequencies of the dipole modes 
coincide therefore with the bare frequencies Ua of the trap. In experiments this feature is 
useful to measure and calibrate the trap frequencies. 

The first order terms rjl^J (t) together with the displacement amplitudes (b (t) are connected 
to density oscillations with s-wave and rf-wave symmetry. Second order displacement am- 
plitudes hke rj^^fj^ (t) are connected to the octupolar collective density excitations Sn^^bc (^) • 
These and even higher order modes we shall consider elsewhere. 

It follows directly from fl7^ that only certain linear combinations of the first order dis- 
placement amplitudes r^^),'' {t) together with the dilatation amplitudes Cb {t) couple to s-wave 
and d-wave symmetry density oscillation amplitudes 5n^^'> (t) and Sn^^ij (t), while the homo- 
geneous zero order displacement amplitudes i]^^ (t) couple to the dipole modes: 



(81) 




6n (r, t) = 2no ■ 



a<b XaXi, 



Pabity-an 



+^Px{t)rx + ^Py{t)ry + j2pz{t)rz 

''x '^^y ^^z 
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, where 



a,be {x,y,z} 

Pab[t) = dab — T- 



(82) 




L\t)\ - (1 - Sab) XaXb \^vll\t) + 



P00{t) =^Paa{t) 



a 



Pait) 



Va 



Particle number conservation implies the solenoidal constraint 



(83) 



, so that there follows immediately 



P00{t) = ^^Paait) 



(At) , Cy{t) , Cz{t) 



(84) 



a 




In general, a BEC cloud may get excited by a combination of actions, involving transla- 
tions of the trap minimum, rotations of the trap axes, or changes of the curvature of the trap. 
The homogeneous displacement amplitudes r]ji^ (t) correspond to infinitesimal translations 
of the position of the center of the BEC cloud. Thus, a density oscillation with a dipolar 
p-wave symmetry proportional to Pa{t) can be excited by a translation of the minimum 
of the trap. Being mainly interested in the effect of interactions on the collective modes, 
however, we set in the following without loss of generality r]ji^ (t) = 0. With regard to the 
first order off diagonal displacement amplitudes rj^ij (t) we easily identify anti-symmetric 
displacement amplitudes r]^^fj (t) = —r]'^'^J (t) as infinitesimal rotations around a rotation axis 
perpendicular to the rar;,-plane, while symmetric off diagonal amplitudes rj^^fj (t) = r]\^^ (t) 
correspond to transverse shear. So, in the geometry under consideration density oscillations 
with a quadrupolar rfa^-symmetry proportional to Pab{t) can be excited by rotations or by 
transversal shear of the trap. Density oscillations displaying an isotropic s— wave symmetry 
proportional to poo{t), and also showing a quadrupolar (i^-wave or (i^2_j^2-wave symmetry 
proportional to certain linear combinations of the diagonal amphtudespaa(^), can be excited 
by a sudden change of the curvature of the trap. 

According to (1791) the time derivative fi-^ds (r, t) of a phase fluctuation associated with 
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such a density fluctuation 5n (r, t) is given by 



(85) 



. (9 ^ , , 



[I - ££,) Utf (r) Tc 



(m ■ Vr)' I rf^r'^n^F (r') < 



'1 - Ed) Utf (r) 



-3£^ (m-Vr)'i; L dV^nTi.(r') 



Jtf Ir-r'l' 

The three-dimensional integrals over the Thomas- Fermi ellipsoid J^tf we express now as 
one- dimensional integrals using Chandrasekhar's integrals (123]) and (12^ : 



47r r - r 



-riTF (r') = no$i (r) 



(86) 



and also 



(87) 



1 

An 



d r 



r — r 



° 4 Sr. ^ ^ ^ 



The crucial trick to proof this representation for the first moment of the Thomas-Fermi 
density profile htf (r) is to use the identity 



-no- 



\ld_ 
4 dr. 



1 



' X y ' z 
A2 A2 A2 

'a; 'j/ '2; 



(88) 



One finds then upon partial integration a surface integral over the boundary ff^TF , and a 
volume integral over the Thomas- Fermi domain 3tf- However, the surface integral vanishes 
identically taking into account that utf (r') = for r' G (9Dtf- So only the volume integral 
contributes, confirming the result ( IHTI) . 

The wave equation (1711) for the density fluctuations we rewrite now 



1 



nTF (r) Vlnps (r, t) 
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(89) 



To calculate the term V'^H-^Ss (r,t) we need 



/ = 1,2,... 



(90) 



Straightforward calculations lead to 



(1 - en) ^Pxx (t) + ^pyy {t) + ^p^^ {t) 



-4:gsno < 



+3eD 



cos^(i?o) 



^p.. it) + ^H^p.. it) + 



sin(2T9o) 
2Ax Az 



Pxz it) 



+ 



2 



A2 



J_ 

Xi 



I 3e_D / cos2(,9o) I sin2(,9o) 



Poo it) 



Likewise we obtain 



a,b,ce {x,y,z} 



driTFiTc) d d 

-h—5s ir,t) 



^ dr„ dtn dt 



(1 - en) Ea,c ^ ^V^aJ it) + it) 



1 ^(1) 



3p. V 9 
Z^a,b,c 4A2 dra 



vi'J it) X 



2 92 



c drijdrc 



>■ Vr)'$2 (r) 



-(1-^d)E..¥^(1 + 25, 



ab TT 



a 

dXb 



(m- Vr)'$i(r) 



(91) 



(92) 



A glance at Chandrasekhar's representation fl26|) for the potential functions (r) of inho- 
mogenous ellipsoids reveals, that for a point r inside the ellipsoid 3rpF the potential function 



for / = 1 is a 4-th order multinomial in the variables |l, 
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nrt^ lY'^ 



} with coefficients propor- 



tional to the index integrals la and lab , and for / = 2 it is a 6-th order multinomial with 
coefficients proportional to the index integrals la , lab and labc (see appendix [A]) . The task 
to calculate the term ( 192|) is therefore reduced to calculate linear combinations of derivatives 
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of certain multinomials in the variables {^l,rl,ry,r'^^ 



a,b,ce {x,y,z} 



b,c 



drbdvc 



(93) 



and 



F,.,. (t) rl + Fyy (t) rl + F,, (t) 



a,b,ce {x,y,z} 



l + A)(m.Vr)^..(r) 



(94) 



Gx. (t) r2 + G,, (t) r2 + G.. (t) 
+Gxj; (t) r^ry + Gy^ (t) r^r^ + (t) r^r^ 



Gxy — — ^yz 



In terms of the triple index integrals (see appendix [XI) 



labc 

(AxjAjyjA^) — A^AjyAz / 



(95) 



^(A2 + (A2 + u) (A2 + «) + (\ + U) (A? + ^) 



the coefficients Fab{t) are determined as linear combinations of the displacement fluctuation 
amplitudes "riabi^) ' ^^"^ coefficients Gab{t) are determined as linear combination of the 
dilatation fluctuation amplitudes Ca{t)- To evaluate the gradient terms in the wave equation 
fl5^ . however, only the differences Gab{t) — Fah{t) are needed, which can be represented as 
linear combinations of the fluctuation amplitudes pab {t) deflned in (152]) . Explicit expressions 
for Gab{t) — Fab{t) in terms of the triple index integrals labc are presented in the appendix 

o 
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Altogether we find for tlie gradient part of (189|) the following second order multinomial 



in the variables |l, 



Tx 



a e {x,y,z} 

-E 



(96) 



nr 



-Woo 



a<b 



ran 



For the coefficient proportional to unity we find 



:i - en) 



-^2 Pxx (t) ~l" y2 Pyy (0 ~^ \2 Pzz if) 



Woo (t) = < 



+3e 



D 



^^P.. it) + ^P.. it) + %g^p.. (t) 



i-sp Mil 



+ tV + t^ + 



A^ 



3eo cos^(i?o) , sin^(i?o) 



A? 



+ 



Poo (t) 



2 

The coefficients of the diagonal terms ^ for a G {x, y, z} in (l96l) are 



(97) 



(98) 



W^aa (t) = Woo (t) + (1 - £d) 4 PP^a (t) + PoO (t)] + ^SdXI [Gaa{t) - F,,(t)] 

while the coefficients of the off diagonal terms for a, G {x, z\ and a < 6 are 



(99) 



Wah (0 = (1 - ( ^ + ^ j Pafe (i!:) + ^e^oAaAb [Ga6(t) " ^^afo (^)] 



It follows directly from what has been said that the right hand side of fl89p represents a 
second order quadratic form: 



iSpoo(t) + ^u'oo(t) 



= < 



„ I 1 9M i 2nog, 



[|.Paa(t) + i^P00(t) + ^«;aa(t) 

+ E, 



7" 



ftpa.(t) + ^^a6(t) 



AaA^ 



(100) 
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Equating the coefficients of the hnearly independent basis functions 1,^-,-^^ for a,b G 
{x, y, z} to zero leads to a set of seven coupled ordinary differential equations for the sought 
fluctuation amplitudes Pabif)- 

As a matter of fact, the equation for the variable poo{t) is obsolete, because the solenoidal 
constraint (Ell) implies 

P00(0 = Pxxit) + pyy{t) + P^^()f:) (101) 



This is consistent because certain identities obeyed by the triple index integrals labc imply 
the following sum rule (see appendix iBj) : 



(102) 



Indeed, adding the differential equations for the diagonal fluctuation amplitudes proportional 
to p- leads immediately to 



(103) 



Consequently the derivative term ^pm{t) and the term WQQ{t) in the differential equations 



fllOOl) for the diagonal density fluctuation amplitudes Paa{t) cancel each other. We obtain 
finally the following six differential equations for six fluctuation amplitudes Pabif)'- 



a,be {x,y,z} 



= T^Paa t + — 



I-Ed)^ 3paa {t) + Eb^a Pbb {t) 



(104) 



Q2 

dt 



2nog 



nr 



;i - eo) (4 + ^) Pab it) 



+ |£DAaAfe [Gabii) — Fab (^)] 



To determine the eigenmodes of oscillation we look for a solution of the form 



Pab{t) = Pab (fi) COS {Qt + 6q) 



(105) 
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, where fl is the eigenfrequency of the mode, and pab denotes a component of the asso- 
ciated eigenvector: 



2nog 





r 


c 

'^xx,zz 


'^xx,xz 










Pxx (^) 




Pxx (^) 




^yy,yy 


^yy,zz 


Cyy^XZ 










Pyy (^) 




pyy (^) 


'^zz,xx 


'^zz,yy 


^ zz^zz 


c 

'^zz,xz 










Pzz (^) 




Pzz (fi) 


Cxz,xx 


Cxz,yy 


Cxz,zz 


Cxz,xz 










Pxz (fi) 




Pxz (fi) 














Cyz,yz 


Cyz^xy 




Pyz i^) 




Pyz (fi) 














Cxy,yz 


Cxy,xy _ 




_ pxy (fi) _ 




Pxy i^) 



(106) 



We find it convenient to ehminate the interaction constant using ([3 

Js) ,,2\2 



2n, 



9^ 



^y\ 



(107) 



1 - £b + [cos2 (l?o) + sin^ (4) 4^] 

It follows upon inspection of the coupled differential equations f ll04p for the coefficients Cab,cd 
explicit expressions, which are listed in the appendix O 

The collective modes associated with the 4 x 4-sub matrix in fll06p describe small am- 
plitude oscillations of the density, which are linear combinations of s-wave and quadrupolar 
dx^-y^ , and ^^.^-waves, while the modes associated with the 2 x 2-sub matrix describe 
small amplitude oscillations of the density consisting solely of combinations of quadrupolar 
dyz-Biid dxy-'waves: 

^Pxxi^yi + ^pyyi^yi + ^fizi^yi 

Poo{^) 



Snn (r,t) = 2no 



1 



r 



' z 



' a 



cos {Qt + 



(108) 



By construction there holds 



d^r Snn (r, t) = 



It is instructive to visualize the spatial dependence of the eigenmodes of small amplitude os- 
cillations of the density by plotting the instantaneous boundary of the BEC cloud when only 
a single mode with eigenfrequency Q is excited. This instantaneous boundary is implicitely 
defined as the surface 

nn{r,t) = nTFir) + 5nn{r,t) = (109) 
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Finally, let us discuss which collective modes can be excited by changing the trap poten- 
tial, always keeping the trap strictly harmonic while changing it. It follows directly from 



a,be {x,y,z} 



(110) 



Pab(^) = Sab 




Sudden changes of the trap potential may excite collective density oscillations around the 
quantum degenerate groundstate. For example, a rotation around a trap axis perpendicular 
to the a6-plane, as represented by the anti-symmetric components of the tensor rj^J , or 
changes of the curvature of the trap, as represented by dilatation amplitudes Ca (^), but also 
transversal or longitudinal shear movements of the trap, as represented by the symmetric 
components of the tensor rf^J, can be used to excite the collective modes f lHTj) of the particle 
density of a trapped BEC cloud. A sudden translation of the origin of a harmonic trap, 
on the other hand, only excites the dipole modes with eigenfrequency Qa = ^a- It should 
be noted, that during these collective oscillations of a spin polarized dipolar BEC cloud, as 
described by the density fluctuation f llOSp . the atoms always keep the orientation of their 
magnetic moments strictly along the external polarizing field B. 



B. Pure Scissors Modes and Mixed Monopole- Quadrupole Excitations. 

Consider a harmonic trap where the principal axis t of the trap is aligned parallel to 
the polarizing external field B, i.e. i3t = 0. In this case the off diagonal matrix elements 
Cxy,yz iCyz^xy , C^z.aa and Caa,xz vauish identically for arbitrary strength Sd oi the dipole 
interaction parameter. There follows then a simpler eigenvalue problem determining the 
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eigenmodes of the small amplitude density oscillations: 



i3t = (111) 
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^xx,zz 













Pxx {^) 




Pxx {^) 




Cyy^XX 


^yy,yy 


Cyy^ZZ 













Pyy (^) 




pyy (^) 




Czz,xx 


Czz,yy 


Czz,zz 













Pzz (^) 




Pzz m 













Cxz,xz 










Pxz {^) 




Pxz {^) 
















^yz.yz 







Pyz (^) 




Pyz (^) 



















'^xy,xy 




. Pxy i^) _ 




. pxy (fi) _ 



Three modes with indices a b display a pure quadrupolar d^z-, dyz- and (ia-^^-symmetry. 
Also there exists a mixed symmetry coupling between two basis functions with d-wave sym- 
metry and one basis function with s-wave symmetry. This is reminiscent of the symmetry 
of the discrete group D^h lifting the 5-fold degeneracy of the / = 2 spherical harmonics into 
three one-dimensional manifolds, namely Aig , Big and B2g- , and a two-dimensional Eg- 
manifold. The one- dimensional (trivial) representation of the isotropic basis function with 
s-wave symmetry we refer to as aig. 

In the geometry under consideration the iJ^-manifold is spanned by basis functions with 
dyz- and d^^y-symmetry, while B2g is spanned by a single basis function with d^^-symmetry, 
and Big is spanned by a single basis function with (i2,2_j^2-symmetry. The one-dimensional 
manifold Aig represents a fixed linear combination of basis elements with ^^2- and s-wave 
symmetry. So, the upper 3x3 block in fillip describes a coupling between members of 
the aig , y4ig-and Big- manifolds. For uj^ = ojy ^ oJz there exists a pure -Bi^-mode, and two 
coupled modes with mixed aig- and Ai^-symmetry. 

The eigenfrequencies of the i?2g-and Eg-modes are obtained from the diagonal matrix 
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elements Cxz,xz , Cyz,yz and Cxy,xy taking the limit i9o 0: 



^z J 1 S£) -\- zy 



(112) 



\l\ (1 - Ed) + Iso'-^Iyzz 



1 tD 2 \2J-zy 



3 A? A2- 



A2\ {l-eD) + hD^^h 



2 xi^'^y 

The spatial variation of the associated density fluctuation of these modes is purely two- 
dimensional 

a, 6 G {x, y, z} (113) 

Pa'b' (^afe) = (^aa'^feft' 

(r, t) = 2nopab{^ab)-r^ cos (l^^fe^ + (^f^^J 

In the limit 5n — > it is found that ^ ^ , Then one obtains for a BEC without 
dipole-dipole interactions confined inside a harmonic trap: 

a^b (114) 



lim Qab = \ Lol + Lol 



These are the so called "scissors" modes first predicted by Guery-Odelin and Stringari 12| . 
and then observed in experiment 

, Q. 

In order to specify conditions that enable excitation of the scissors modes (11131) for 
a dipolar BEC cloud confined in a harmonic trap we point out, that the components 
Pa'b' i^ab) = ^aa'hb' of the eigenvectors of the respective modes are connected to the off 
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diagonal displacement amplitudes ry^^ , see f l82|) . by 



a^b (115) 

Pab (^ab) = — AaAb 

For an infinitesimal rotation of the BEC cloud around one of its symmetry axes, say e^fi = 
Ga.oA ebfi , the associated displacement amplitude is anti-symmetric, rf^J = — ^^P- So one 
recognizes immediately that in the highly symmetric case "(^t = a scissors mode with 
amplitude 'pab {^ab) cannot be excited by a rotation around a principal axis of the BEC 
cloud perpendicular to the a6-plane, if the semi-axes Xa and Xb of the BEC cloud in that 
plane are equal, i.e. = A^. However, even then a scissors mode may get excited by 
a sudden tranverse shear movement of the trap as described by a symmetric displacement 
amplitude rj^J = rf^a ■ BEC is confined inside a harmonic trap with tri-axial symmetry, 

one may always excite the scissors modes Vtab by a sudden infinitesimal rotation of the trap 
potential around a symmetry axis perpendicular to the respective a6-plane. 

Let us now discuss the coupled modes corresponding to the 3 x 3-sub block in f illip . 
These are small amplitude oscillations of the density that are linear combinations of the 
three diagonal amplitudes In the limit — > the corresponding eigenfrequencies 

and eigenvectors of the triplet of coupled modes can be obtained solving a cubic equation 
for the frequencies 

(116) 
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(fiW) " 




Pxx 


(fiW) " 






Pyy 


(^](0)) 


Pyy 


(^](0)) 






Pzz 


(fiW) 




Pzz 


(fiW) 



One easily sees, that for a tri-axial trap the eigenmodes of this triplet are mixtures of 
basis functions with isotropic s-wave and quadrupolar dz^-v^me and dx2_y2-wme symmetry, 
respectively. 

When the harmonic trap has a uniaxial (cylindrical) symmetry, Uz ^ ^y = ^x = , 
simple analytic formulas for the eigenfrequencies and eigenmodes of the density oscillations 
of a BEC cloud can be derived from flll6p that apply for ££> = 0. One easily obtains for 
the three eigenfrequencies fl^^2_ 2 , ^f''^ , ^^^2 well known results first derived by Stringari 

X y z 



ab) 



34 



In the appendix [D] we present a detailed discussion of these modes as a function of the 
anisotropy ratio 

(117) 



C. Spherical Harmonic Trap 

For the special case of a spherical harmonic trap, say with trap frequency Ua = oj, setting 
aI^^ = A, we immediately find from f lll6p for a BEC without dipole-dipole interaction (see 
appendix [D]) : 

n = ^2^2,2 = V2u (118) 



6nQ (r, t) = 2no cos {Qt + 6, 



r1 — 



y 



A2 



Srin (r, t) = no cos {Qt + 6n) 



(119) 



' X ~^ y ^ z 

' A2 



^] = = y^oj (120) 



<5?2n (r, t) = no cos (fit + Sq) — - — -| 

So for = 0, a BEC cloud confined inside a harmonic spherical trap may get excited as an 
s-wave breather mode with frequency fi^"* = y/Eu , or as a quintuplet of degenerate modes 
with quadrupolar symmetry and frequency ^^^2_y2 = ^^^^ = ^15 = ^yz = ^xy = V^u, 
namely three scissors modes with d^z-, dy^-, (iajj^-symmetry, and two modes with d^^^yi- and 
(i;22-symmetry. 

Next we take into account the effect of the dipole-dipole interaction. According to (149|) 
for Ed > the groundstate of a spin-polarized dipolar BEC cloud confined in a spherical trap 
with trap frequency Ua = uj displays uni-axial symmetry along the direction of the magnetic 
field B, so that A^. = A^ < A^. Let us check if for eo the modes of a dipolar BEC cloud 
confined in a spherical trap are qualitatively similar to the aforementioned collective modes 
of a BEC cloud without dipole-dipole interaction, ed = 0, for the case of a prolate trap with 
cylindrical sjTiimetry: Uz < ojy = Ux- 
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Indeed, for a spherical trap with trap frequency u we have i9o = , so that all matrix 



elements in (lllip can be expressed in terms of the following expressions: 

1 



It follows then from f illip : 



A 



B 



C 



' 2^ ' yy^ 



2 A2^2;y 



A,, 



X'i ^ 2/22 



1 — ^r, -I- 
1 c- -J- ^Eg -''hT 



yz 



;i2i) 



(122) 



Pa'b' i^ab) — Saa'^bb' 

So, the scissors modes with rfaj^-and d^^-symmetry remain degenerate. 

For 7^ the 3x3 sub block in f illip represents a triplet of coupled modes. For the 
case of a dipolar BEC confined in a spherical trap there follows 



3A A 
A 3A 



B 
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^B ^B 3^C 

A~ A, A, 
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Pxx (^) 




Pyy (^) 




Pro (^) 




(^) 




Pzz i^) 



It is easy to see that the mode with quadrupolar d, 
eigenstate for 7^ 0: 



x^—y^ 



(123) 



symmetry remains an exact 



xy 



(124) 



Pxx iS^x^—y'^ ) 
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Pyy iS^x'^—y'^ ) 




-1 


Pzz iS^x^—y'^ ) 
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So, for 5d 7^ the quadrupolar modes with dx2-y2-and dxy-symmetry remain degenerate for 
the case of a spherical harmonic trap. 

Next we show, that the isotropic breather mode of a dipolar BEC inside a spherical 
trap with frequency u is an eigenstate of the small amplitude density oscillations of the 
BEC cloud, displaying an exact (!) s-wave symmetry for any value of the dipole interaction 
strength 7^ : 



(125) 



Pxx i^s) 
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Pyy (^s) 
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Pzz iP's) 
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If this claim was correct then it should be true that 



AA + B 
{2B + 3C) 



5 
5 



(126) 



Indeed, making use of identities (lATp . flASP obeyed by the triple index integrals lahc we see 
that 



4—1 +37 



yzz 



5/ 



zy 



(127) 



4y4 + S = 4- 



1 — + 9^1) ■ -^lyyz 1 — £d + n^D ' tI-^ 



i- £D + —Xljzy 



+ 



1-60 + —Xljzy 



5(1 -£d) + f£D 



3_ _ f.X^ 



A? 



4 _i^2 -^j/j/^ + 3/ 



yzz 



1 c- L 3£r) "^y r 

i - fcD + —J^J^zy 
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, and also 



zzz ~\~ '^'T^Iyzz — zz (12^ 



Kr Xl^-eD + ■ ^lyzz XI 1-60 + ^EdLzz 

-I (25 + 3C) = 2-f \ + 3-f ' 

A2 5 (1 - sd) + I^D (2§lyzz + 57. 
y \ ^ 

Z L t£) -T 9 \2 J- zy 

5 ^ 



.Xy 1 — Ed + ^EdIzz 

1 _ I 3£D ^7 



The last line follows because the selfconsistency equation ( H6l) implies for the case of a 
spherical trap: 



\2 1 _ <r„ -I- ^ r 
\ _ tD -t- 2 \i^yz 



(129) 



Because for a finite value 5^ > a spin-polarized dipolar BEC cloud confined in a spherical 
harmonic trap with frequency u has the shape of an uniaxial (prolate) ellipsoid orientated 
parallel to B , so that A^. = Ay < A. , we find it remarkable that the isotropic breather mode 
(11251) remains an exact eigenmode of the small amplitude density fiuctuations with s-wave 
symmetry, oscilUating at a constant frequency Vts = y/Eu that is independent on the value 
of the dipole interaction strength for — | < ed < 1- 

Knowledge of two eigenvalues is sufficient to determine the third one from the trace of 
the coefficient matrix in (I123p : 



+ nl,_y, + nl, = qu' (^a + (i30) 



This leads for the eigenfrequency and the eigenvector of the density oscillations with a 
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predominant (i^2-symmetry to the result: 

J 



= {4:A + 3^C -5] (131) 



5-4A 
XI 2B 



5-4A 
A2 25 



Pyy (^2 
P22 (^2 



It follows from what has been said that the degeneracy of the small amplitude collective 
modes of a dipolar BEC cloud confined in a spherical harmonic trap is only partially lifted 
for Ed ^ Q. For a spherical trap the modes with (i^.2_y2-and with (ia-^^-symmetry, and also 
the modes with dyz- and daj^-symmetry remain degenerate, irrespective of the value of the 
dipole interaction er,. In FigJHlwe plot the collective mode frequencies VLg ,^22 , ^xy and VL^z 
vs. the interaction strength parameter ed- For small \ed\ the splitting of the quadrupolar 
modes ^2^2 , Vt^y and Vtxz is weak. Most remarkably, the breather mode Vts displays for 
— 1 an exact s-wave symmetry, the eigenfrequency assuming a constant value 

VLs = \f^uj, even though for s/j 7^ the shape of the groundstate is not isotropic. 

The following reason can be given for the breather mode frequency of a dipolar BEC 
being independent on the dipole interaction strength Ed for an isotropic harmonic trap. 
The microscopic Hamiltonian of a dipolar interacting gas cloud consisting of atoms is 



H = H, 



kin T J^pot 



H„ot + H,. 



int 



(132) 



2m* 

"=1 aG{x,y,z} 



N 



*, ,2 



E ^^^^ 



pot 



n) 



n=l 



aE{x,y,z} 



int 



n,n'=l 
n'j^n 



The breather mode (or monopole mode) of small amplitude collective density oscillations 
of such an atom cloud may get excited by a sudden change of the curvature of the trap 
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potential, say by changing the trap frequency w — )■ a; + 5u. The associated excitation 
operator is 

N 

^ ^ " " (133) 



6V = m*u6uJ2 

n=l a£{x,y,z} 

It is important to reahze, that the interaction potential U (r*^"\ r*^"')) for the spin-polarized 
dipolar BEC in transforms under a scaling transformation r — Ar like a homogeneous 
function with scaling degree —3: 

(134) 



Together with Newton's law of action and reaction fl78|) this implies 



dV.H 



2ffu6u) {2Hpot-2H, 



' kin 



3K 



int 



(135) 



If ^'o denotes the groundstate and Eq the groundstate energy of the system under consider- 



ation, there holds 



0= (^0, 

^0, 



6V,H 
6V,H 



Eq — Eq 



dV,H 



(136) 



H 



^0 



Inserting the double commutator fll35p it is found, that the full interaction energy of a 
dipolar interacting BEC in the groundstate is proportional to a difference of kinetic and 
potential energy only: 



int 



(137) 



*o 3 

It should be emphasized, that if in (jl]) the scaling degree of the long ranged interaction ([6]) 
under r — t^A r was different from the scaling degree —3 of the short ranged s-wave contact 



would not apply! 



interaction (|5]), the derived virial identity ([13 

Next we employ a well known sum rule 10| , ll| providing an upper bound for the low- 
lying excitation energies Ei — Eq that can be excited by a hermitean perturbation operator 
6V : 



{E, -EqY< 



For the operator 5V exciting the breather mode, see (I133p . it is found 



(^0, 


6V,H 
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6V,H 
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H,6V 



8h' {Suf Hpot 



(139) 
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and also 



(140) 



5V,H 



6V,H 



H 



From what has been said there follows now for the frequency Qs of the breather mode an 
upper bound: 



*0 



2 ( Hpot ) ^ + 2 ( Hkin ) _ + I ( Hint 



*0 



pot 



(141) 
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kin 
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For the optimized groundstate ([T]) of a BEC, as constructed from a solution to the Gross- 
Pitaevskii equation ([2]), the ratio of kinetic to potential energy scales like 



\ / *o 



pot 



o N--^ 



(142) 



So, in the Thomas-Fermi approximation the derived upper bound for the breather mode 
frequency is indeed independent on the strength of the dipole-dipole interaction parameter 
Eo- The fact, that this upper bound actually coincides with the previously derived result 
Vlg = a/Sco' , which was obtained solving the eigenvalue problem fll23p for the small amplitude 
collective modes of density oscillations, suggests that the spectral weight of the mode is 
indeed exhausted by the specified excitation operator 6V (1133^ of the monopole mode. 

It is instructive to visualize the spatial variation of the associated density eigenmodes 
6nQ (r,t) by plotting the instantaneous surface of the BEC cloud as defined by fllOQp . In 
Fig. 7 and in Fig. 8 these eigenmodes are plotted at stroboscopic times t = , t = ^ and 
corresponding to maximal, zero and minimal deviation from the boundary OBitf 



of the groundstate cloud Dtf , respectively. The plots shown are based on a selfconsistent 
calculation of the groundstate cloud for a dipole interaction strength parameter £_d = 0.7 
, assuming that the BEC cloud is confined inside a spherical harmonic trap with trap fre- 
quency u. The amplitudes of the respective eigenmodes 6n^ (r, t) of the density fluctuation 
have been scaled by a suitable factor for each mode separately to make the typical shapes 
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FIG. 6. (Color online) Eigenfrequencies Qg , > ^xz = ^yz and ^l^y = ri^2_y2 of small amplitude 
density oscillations of BEC cloud vs. dipole interaction strength e^i for a spherical harmonic trap 
with trap frequeny uj. 

better visible. The s-wave breather mode is clearly distinguished in its appearance from the 
three characteristic scissors modes with their d^z , dyz and d^y-^SNe symmetry, and also the 
(i2,2_j^2-wave and (i^-wave quadrupolar modes. 
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FIG. 7. (Color online) Visualization of density fluctuations nn(r, t) = nTri^) + ^n^i^^, t) of scissors 
modes for dipolar BEC cloud confined inside a spherical trap for a dipole interaction strength 

ED = 0.7. 
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t-n = o 




n^2_y2 = 1.30022c.; 
t-n = n/2 




t-n = o 




= 1.4011 w 
t-n = 7v/2 



t • 17 = TT 




FIG. 8. (Color online) Visualization of density fluctuations ^^(r, t) = nTpi^^) + '^^.^(r, t) of dipolar 
BEC cloud confined inside a spherical trap for a dipole interaction strength £d = 0.7. First row 
isotropic breather mode ilg ; second row quadrupolar mode Q^2_y2 ; third row quadrupolar mode 

r2~2 . 
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D. Spectrum of Low-Lying Excitations for the Case 'dx = 0. 



We now discuss the collective density oscillations of a dipolar BEC cloud confined in a 
tri-axial harmonic trap in the highly symmetric case, when the principal axis ^ of the 
trap is orientated colinear to the spin polarizing magnetic field B, so that dx = In Fig. 9 
and Fig. 10 the collective mode frequencies VL corresponding to the solution of the eigenvalue 
problem f illip are plotted vs. the dipole interaction strength ed- Shown are three scissors 
modes with d^z , dy^ and d^^-wave symmetry, and three hybridized modes combined from 
basis elements with s-wave, dz^-wave and dx2_y2-vfaMe symmetry. The anisotropy ratio chosen 
is Ux '■ ujy '■ ujz = 712 : 128 : 942 in Fig. 9, and in reversed order UJx '■ UJy '. bJz = 942 : 128 : 712 
in Fig. 10, respectively. 

: ojy : oj, - 712 : 128 : 942 , &t = 0° 
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FIG. 9. (Color online) Eigenfrequencies of small amplitude collective modes combining isotropic 
s-wave and quadrupolar d-wave basis elements vs. dipole interaction strength ed for dipolar BEC 
cloud confined in a harmonic trap with anisotropy ratio cOx ■ ujy : uJz = 712 : 128 : 942 in the highly 
symmetric case = 0. Displayed are three scissors modes (blue lines) and three hybridized modes 
that are combinations of s-wave, dx2_y2-a.nd d22-basis elements (red line, black line and green line). 

There exists fair a gree ment between our exact analytical results and the numerical re- 



sults obtained in Ref. 



21| . which are based on the method of solving Newton equations of 



motion for time dependent Thomas- Fermi radii. As is evident from (1 11 01) . small amplitude 
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cj^ : cjy : cj. = 942 : 128 : 712 , (^r = 0° 
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FIG. 10. (Color online) Dependence of eigenfrequencies of small amplitude collective modes com- 
bining isotropic s-wave and quadrupolar d-wave basis elements vs. dipole interaction strength ed 
for dipolar BEC cloud confined in a harmonic trap with reversed anisotropy ratio uJx : ujy : uJz = 
942 : 128 : 712 in the highly symmetric case tDt = 0. Displayed are three scissors modes (blue 
lines) and three hybridized modes that are combinations of s-wave, dx2_y2-and (i^2-basis elements 
(red line, black line and green line). 

fluctuations of the Thomas- Fermi radii are described in our approach by the dilatation am- 
plitudes Ca{t)- However, in the highly symmetric case -i^T = 0, these dilatation amplitudes 
only couple to the diagonal basis elements of the tensor pab '■ 



An 



(143) 



In the highly symmetric case = no coupling of the dilation amplitudes Ca (^) to the off 
diagonal elements a 7^ 6 of the tensor pab exists, as is evident from (11151) . To ease comparison 



of our results with the results presented in Ref.( 



we also plot in Fig. 11 and Fig. 12 the 



relative change of the collective mode frequencies vs. Ed for the three hybridized 

modes displayed in Fig. 9 and Fig. 10 that couple via the dilatation amplitudes Ca to the time 
dependent Thomas-Fermi radii. 

It should be pointed out that a purely diagonal shear movement of the dipolar BEC cloud 
at constant Thomas-Fermi radii, Ca = 0, as described by the diagonal elements rfaJ of the 
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cj^ : cjy : = 712 : 128 : 942 , &t = Q° 
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FIG. 11. (Color online) Relative change — ^ — vs. ed for the three hybridized collective modes as 
displayed in Fig. [9] for the highly symmetric case i9r = 0. 

oj^ : (jjy : cj^ - 942 : 128 : 712 , (^r = 0° 

i-n(°Vn 




FIG. 12. (Color online) Relative change — ^ — vs. ed for the three hybridized collective modes as 
displayed in Fig. [10] for the highly symmetric case '&t = 0. 

tensor rf^^ spanning the (solenoidal) displacement vectorfield (IHOj) . may also excite these 
hybridized modes coupling to paa- This degeneracy is a special property of any quantum 
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degenerate BEC groundstate with an ellipsoidal shaped density profile. 
E. Spectrum of Low-Lying Excitations for the Case -dr / 0. 

Sudden changes of the trap potential may excite various collective modes of a dipolar BEC 
cloud. If the polarizing external magnetic field B is not in alignment with the principal axis 
e^^T of the trap, so that ^ includes a finite angle 'Ot ^ ^ with B in the xz-plane (see 
Fig{T]), the s-wave and d-wave symmetry parts of the collective density oscillations combine 
to a quadruplet and a doublet of modes. It is found from (11061) . that the modes with mixed 
dx^-y^ , dz^, dxz- and s-wave symmetry, consisting of a linear combination of the three 
diagonal amplitudes p^^ , Pyy , Pzz and one off-diagonal amplitude p^z, combine together 
to a quadruplet (see Fig. 13), and the modes with mixed d^y- and liy^-symmetry combine 
together to a doublet of scissors modes (see Fig. 14). From (11151) it is evident, that for 

7^ an infinitesimal rotation around the principal axis e^^T of a harmonic tri-axial trap 
may then excite via its coupling to the paj^-components of the eigenvectors all four modes 
of the mentioned quadruplet of small amplitude oscillations of the density simultaneously. 

Likewise, a rotation around the principal axis e^, (or e^. ^ ) of the tri-axial harmonic 
trap may excite via the coupling to the off diagonal amplitudes 'pyz and 'pxy the mentioned 
doublet of scissors modes simultaneously. Alternatively, these scissors modes can also be 
excited by transversal shear movements of the anisotropic harmonic trap, thus creating an 
excitation of the BEC cloud that may be described by a (solenoidal) displacement vectorfield 
( IHOl) that is spanned by the symmetric off diagonal elements of the tensor rf^fj . 

A sudden change of the curvature of the trap potential, as described by the dilatation am- 
plitudes (a in (I143P .excites in the geometry under consideration the modes of the quadruplet, 
but never the scissors modes of the doublet with mixed d^y- and rf^^-symmetry. 

The results displayed in Fig. 13 and Fig. 14 reveal, that a tri-axial harmonic trap with trap 
frequencies = Ui , Uy = U2 , ujz = (jJs , say ui > U2 > ujs , shows a characteristic shift of 
the eigenfrequencies of these quadruplet- and doublet-collective modes compared to a trap 
with reversed trap frequencies, i.e. a harmonic trap with Ux = uj^ , Uy = i02 , uiz = 

From measurements of these characteristic shifts of the collective mode frequencies of 
the quadruplet- and doublet-collective modes of a dipolar BEC cloud for two such mutually 
reciprocal tri-axial traps the strength of the interaction parameter Sd could be determined 
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FIG. 13. (Color online) Dependence on dipole interaction strength e_D of eigenfrequencies of 
small amplitude density oscillations corresponding to 4 x 4-block in (jl06p when the BEC cloud is 
confined in tri-axial harmonic anisotropic trap for cj^,' : f^j/ : f^z = 6 : 3 : 2 (red line) and 
= 2:3:6 (black line), choosing a trap orientation angle tDt = 5.7°. All frequencies normalized to 
geometric mean uj = {ujxUjyUJz)^ . 



accurately 2l(] . Knowing the mass m* and the magnetic dipole moment |(M)| of a single 
atom one then obtains immediately from ([8]) the isotropic s-wave scattering length of the 



atoms 



/"o|(M)r 



127rri2 



(144) 



The experiment suggested here consists in preparing a quantum degenerate spin polarized 
dipolar BEC cloud confined in a harmonic trap with tri-axial symmetry, so that the principal 
axis 5" of the trap is first orientated colinear to the spin polarizing magnetic field B, i.e. at 
the beginning of the experiment = ^ = (see Fig{T]). Then, say at time t = 0, the trap 
orientation angle is changed suddenly to a new value, by making a rotation around the 
principal axis Gy^T of the trap by a constant small rotation angle, say = 5.7° , the value 
chosen in Fig. 13 and Fig. 14. A dipolar BEC cloud excited in this manner will then oscillate 
not around the old cloud orientation angle -i^o = , but around a new cloud orientation angle 
'Oq {'dx), which is via the self consistency equations (H5l) . (j46l) . (j47|l not only dependent on the 
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CD., : cOy : CO, = 6 : 3 : 2 & 2 : 3 : 6 , &t ^ 5.7° 
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FIG. 14. (Color online) Dependence on dipole interaction strength e_D of eigenfrequencies of 
small amplitude density oscillations corresponding to 2 x 2-block in (jl06p when the BEC cloud is 
confined in tri-axial harmonic anisotropic trap for cj^,' : f^j/ : f^z = 6 : 3 : 2 (red line) and 
= 2:3:6 (black line), choosing a trap orientation angle tDt = 5.7°. All frequencies normalized to 
geometric mean uj = {ujxUjyUJz)^ . 

strength of the dipole interaction parameter ed , but also on the chosen trap orientation 
angle ^9t- The principal axis o of the new equilibrium BEC cloud confined in a harmonic 
trap with trap orientation angle 1?^ then includes with the fixed magnetic field B a finite 
angle "^o, that is smaller or larger than 'Ot, depending on the anisotropy ratio of the trap 
(see Fig.2 ). It follows from what has been said that the eigenfrequencies Q of the collective 
modes of the density fluctuations, that can be excited in this manner, are functions of Ed 
and the trap orientation angle -^t- 

In Fig{T5] and Fig{T6] the dependence of the collective mode frequencies of a dipolar BEC 
cloud on the trap orientation angle t?T is shown for three values of the interaction strength 
parameter Ejj. 
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FIG. 15. (Color online) Dependence on trap orientation angle 'St of eigenfrequencies of small 
amplitude density oscillations corresponding to 4 x 4-block in p06p when the BEC cloud is confined 



in tri-axial harmonic anisotropic trap for uj^ : cOy 



6:3:2. All frequencies normalized to 



geometric mean uo = (w^w. 



The strength of the dipole interaction is = 0.2 (dotted line) , 



Ed = 0.5 (dashed line) , = 0.8 (solid line). 
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FIG. 16. (Color online) Dependence on trap orientation angle 'St of eigenfrequencies of small 
amplitude density oscillations corresponding to 2 x 2-block in (jl06p when the BEC cloud is confined 
in tri-axial harmonic anisotropic trap for uj^ : cOy : uj^ = Q ■ 3 : 2. All frequencies normalized to 



geometric mean uj = {uJxUJ' 



The strength of the dipole interaction is = 0.2 (dotted line) , 



Ed = 0.5 (dashed line) , ed = 0.8 (solid line). 



If a dipolar BEC cloud is confined in a harmonic trap with tri-axial symmetry, so that the 
spin polarizing magnetic field B is orientated in a completely general fashion, i.e. B its not 
orientated parallel to any symmetry plane of the trap, the quantum degenerate groundstate 
of the BEC cloud is then characterized by three Euler angles determining the orientation 
of the ellipsoid Dtf relative to the axes of the trap. A study of such spin polarized dipolar 
BEC clouds we shall present in a separate publication {22!, together with a discussion of 
the octupolar modes of density oscillations, which can be described by fluctuation ampli- 
tudes rj^^l^ (t) associated with a solenoidal vectorfield with a quadratic spatial variation, i.e. 
Va (r, t) = J2b,.c Vale (^) nrc- 



IV. CONCLUSIONS 



We have studied the groundstate and the low-lying collective modes of a dipolar Bose- 
Einstein Condensate for the case that the external magnetic fleld is not necessarily oriented 
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parallel to one of the principal axes of the harmonic anisotropic trap. In particular, we have 
determined the eigenfrequencies of six low-lying collective modes that combine respectivly, 
to a quadruplet and doublet of atom density oscillations with mixed s- and d-wave symmetry, 
and obtained analytical expressions for them. We have found the following results: the mode 
frequencies depend on the dipole interaction parameter in a characteristic way that could 
be used to measure the s-wave scattering length of the atoms accurately. In the special case 
that the harmonic trap is spherical we find the remarkable result that the eigenfrequency of 
the isotropic breather mode does not depend on the dipole interaction strength, even though 
the shape of the condensate does. Thus, this mode could be used as a reference frequency 
for the other collective modes that depend on the dipole interaction strength. A rigorous 
sum rule argument shows that this feature of the breather mode is a consequence of the 
scaling property ( I134p of the interaction potential in a dipolar BEC, and the Thomas-Fermi 
approximation . 
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Appendix A: Index Integrals 



Consider the index integrals 



a,&,c e {x,y,z} 

r 

Jo 
Jo 

r 

{^xi \i ^z) — ^x\^z / 
JO 



(Al) 



du 



du 1 

^{Xl + u) {XI + u) {XI + u) + + 
du 1 

^Xl + u) {XI + u) {XI + u) + ^) + (Ac' + 



These integrals are symmetric under permutations of the indices a,b,c & {x, y, z}. It is also 
evident that 



la- h — - - A^) lah 

lac ~ he = ~ {^a ~ ^b) hbc 

Also index integrals lab and labc are connected by a derivative operation: 

1 d \ 

Iab — {i + '^Sac + S^fec) A 

c^abc 



Xr dXr 



Let us note the identity 



du 



^{Xl + u) {XI + u) {Xl + u)_ 



E 



^(A^ + U) {XI + U) {XI + U) aeix,y,z} ^« + " 



Likewise 



du 



_^{Xl + u) {XI + u) {Xl + u)^l + ''_ 



^{Xl + u) {XI + u) {XI + u) + ^ 



1 



Xl + u ' ^ Xl + u 

" ae{x,y,z} " 



(A2) 



(A3) 



(A4) 



(A5) 
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and 



(A6) 



-2- 



du 



1 1 



2 2 

+ 



+ 5Z A2 + 

ae{x,j/,z} " 



Upon integration with respect to the variable u from to oo there follow now several useful 
identities: 



2= E ^« 

a£{x,?/,z} 

2 



(A7) 



ae{x,j/,2:} 



bca 



Using these relations various useful algebraic connections between the integral , lab , also 



between lab and labc become evident 



15|: 



a,b,ce {x,y,z} 
a ^ b c 



(A8) 



3/aaA^ + labXl + lacX^ — 'ila 



^laaaXa ~\~ ^aabXfy + laacX^. — 5/q, 



3/aabA^ + SlabbXl + IabcX\ — 5Iab 



Making the substitution 



u-^ u = -u' 



(A9) 



55 



we obtain useful scaling relations 



la (Ax, Xy, A.) = Ia[j, ^, ^) = la (AlO) 

1 ^ /Ax A,, A2 \ 1 



lab (Ax, Ay, A^) - j^Iab y-^, = ]y2'^'=''' 

7- (\ \ \ ^ — J_7' = J_T 

^abc l^x, '^j/, ^zj — J^J^abc I "^5 ~K I ~ ]\4 '"^ 



In our calculations we find it convenient to choose A = A^ > 0. 

The task to calculate a double index integrals lab can be reduced to calculating simpler 
single index integrals la- This is enabled by using 



= ^ (All) 



, an immediate consequence of flA7p . Provided Aa 7^ A^ , the integrals J^q can be reduced to 
calculating the simpler integrals and /„ using the identity: 



ae{x,y} (A12) 

^ za ^2 



1 _ ±<k 



So for Ao ^ \z all double index integrals Iza can be reduced to single index integrals /„. 
Carlson [l6| has provided an elegant efficient algorithm based on the well known method 
of the arithmetic-geometric mean to calculate the single index integral la directly, which 
method we highly recommend because of its accuracy and speed |17i] . 

For \a = Xz the right hand side becomes formally undefined. However, in this case we 
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may calculate the integrals Izx and I^y in closed form: 

(A13) 

T -J T Ax- Ay du 1 

lim Iza — 1™ 



Aa ^ Az Aa ^Az /X^ 



lim /^^ = ^ 

Aa? ^ Az 



lim Izy = 



A.7o 

K f°° du 1 
Az 7o 



In the isotropic case A^ = A^; = A^^ : 

- 2 

lim^ lim^ J^^ = lim lim I^x = /(I) = - (A14) 

Ay ^Az Ax ^Az 5 
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Appendix B: Sum Rule 



The sum rule 



follows directly from the defining equations 
integrals: 



J2^aa{t)=5Wooit) (Bl) 
a 

and the properties of the triple index 



(B2) 



3f 
2^D 



COS^ {'do) (34^^ + Ir,y^ + 3Ia:z 

+3 sin^ (i9o) (S^^^x + Lxy + h 

COS^ {'do) {I^y^ + 3Jyy^ + 3/y 



+ Sin^ (-i^o) {"ilxxy + ^Ixyy + ^: 

3 cos^ {'do) (4^^ + Jy^^ + 5/^ 



Pxx{t) 
Pzz{t) 



+ sin^ (^9o) (3/x.x2 + Ixyz + 3^22 
+ sin (2i?o) A^A^. (34^^ + /^^^^^ + 34^^) ^^^(t) 



+5 (1 - Ed) ^Pxx{t) + ^Pyy{t) + ^fzz{t) 

^ l-gp J_ _^ J_ _^ J_A _^ 9£D /' cos^(i?o) _^ siv?{A)) 
2 \^A| y 2 \^ Al Al 



Poo(^) 



There holds the following identity for the index integrals labc- 

a,be {x,y,z} 



(B3) 



"^laab ~l~ '^labb ~\~ ^abx ~\~ ^aby ~l~ ^abz \2 \ 2 



2 



Taking into account that the index integrals labc are invariant under permutations of their 
indices a,b,c & {x, y, z} it follows for the linear combinations encountered in (1B2|) : 

_ 2 

'^-^xxz ~l~ ^xyz ~l~ ^-^xzz '^^xxz ~l~ '^^xzz ~l~ ^xzx ~l~ ^xzy ~l~ -^xzz ^2 ^2 

2 2 

5 1 XXX ~l~ I XXV ~\~ IxXZ '^IxXX ~l~ '^IxXX ~l~ IxXX ~\~ I XXV ~l~ -^x 



'xxa; I ^xxy 1 -^xxz -^-^xxa; 1 ^-^xxx 1 -^xxx 1 ^xxy 1 -^xxz \2 \ 2 \4 

_ _ 2 

~l~ ""^lyyz ~l~ "^^yzz '^^yyz ~l~ '^■^yzz ~l~ A/zx ~l~ A/zy ~l~ -^1^22 , 2 \ 2 
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xyy I -^xyz 



^■^xxy ~l~ '^■^xvv ~l~ -^xvx ~l~ ^xv'u ~l~ 



xyy I ^xyx i ^xyy i 



'^^zzz ~\~ '^J-zzz ~\~ Izzx ~\~ Izzy ~\~ I, 



Making use of these identities, and taking into account f llOip . we in deed see that 



a&{x,y,z} 



3eD 



+ 
+ 



cos^ (^9o) Jj + 3sin2 (4) p,,(t) 

cos^ (t9o) ^ + sin^ (t^o) ^ Pyyit) 
3cos2 (^?o)^ + sin2 

+ sin(27?o)^ Px.(t) 



+ 



l^Af + Aj + 

sin^(i9o) 



1 A I 9eD cos^()?o) sin^(i3 



Poo(i!:) 



3£ 



D 



+ 



cos^ (t9o) T2 + sin^ (?9o) 4 [p^^(t) + Pyy{t) + p^^(t)] 

+ 5 (1 - Ed) ^Pxxit) + ^ Pyyit) + ^pzz{t) 

Pooit) 



+3eD ■ (3 + 2) 



"f \2" -r 



9£d f cos^(i9o) _|_ sin^(i9o) 
2 \ A? A?. 



sin-'(i?o) 
A~ 



cos (i9o) 



Pzz{t) + 



sin(2i9Q) 

2At; Az 



5- < 



(1 - Sd) ^Pxx {t) + ^pyy [t) + ^pzz {t) 



-3e 



D 



cos-^(i9o) 



-yT^Pzz [t] + "°^2 ° Pxx {t) + 



sin(2'j^o) 

2A'r As 



Pxz {t) 



A^ Xy A| y 2 \ A^ 



Poo (t) 



5woo (t) 
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Appendix C: Explicit Expressions for Gab{t) — Fab{t) and Matrix Elements Cabled 



We collect here explicit expressions for the linear combinations Gabif) — Fabif) that occur 



m 



Gxx{t) — Fxx{t) — — < 



3X1 [cos^ (t?o) Ixxz + 5 sin^ (t^q) Ixxx] pxx (t) 

[COS^ {'do) Ixyz + 3 Sin^ {do) hxy] Pyy {t) 



+3X1 [COS^ (7?o) Ixzz + Sin^ (?9o) Ixxz] Pzz (t) 

+3 sin (2t?o) hxzXxXzPxz (t) 

Xl [COS^ {'do) Ixyz + 3 Sin^ (t^o) Ixxy] Pxx (t) 

r lA^F (A- — ] ^^^y ^yy' + ^^^^ 

yy\ ) yy\ ) \2 i r 

+Xl[3cosH^o)Iyzz + sm^^o) Ixyz pzz {t) 

+ sin (2^9o) Ixyz ■ XxXzPxz (t) 

3X1 [COS^ {^o) Ixzz + Sin^ {^o) Ixxz] Pxx (t) 
^ ^ 1 I +Xl[3cOS^l^o)Iyzz + Sm^'do)Ixyz]Pyyit) 

+3X1 [5 cos2 (^o) + sin^ {^,) p,, (t) 
+3 sin {2^o) IxzzXxXzpxz (t) 



and also in (p9!) : 



Gxz{t) — Fxz{t) — 



3 + If) ^^^^ [COS^ (t9o) + Sin^ (^9o) Pxz (t) 

3 ( 1 + T§ ) IxxzPxx (i) + ( TT + ) Ixyzpyy (i) 

-sin(2^9o) ^ 2 \ ^ ^ 

+3{^^+l)lxzzpzz {t) 



(CI) 



Gyz{t)-Fy,{t)=[^ + ^ 

\ ^z 



sin {2do) XxXylxyzPxy {t) 
+ XyX^ [3 COS^ {do) lyzz + Sin^ {'Qq) Ixyz] Pyz (t) 



XxXy [COS^ {'do) Ixyz + 3 Siv? {'do) Ixxy] Pxy {t) 
+ XyXz sin {2'do) hyzPyz {t) 

The matrix elements Cab,cd occurring in the eigenvalue problem fll06p are explicitely given 
by 
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2no^(") _ 2 ^ (1 - ^i?) + i^o • f [cos^ (??o) Ixxz + 5 sin' (i?o) hxx] 



'Cxx,xx — ^« »2 Q >,2 



l-eD + ^t% [cos2 {do) hy + sin' (^?o) 



m* ""'"^ 1 _ + I [cos2 (^o) + sin' (^q) 4,] 



^* ' 1 - + ^ ^ [cos2 (,?o) Izy + sin' (l?o) 4,] 



2no^ _ 2 |£Dsin(2i?o)^|/..^ 



^ - £^ + [cos2 (l?o) /.J, + Sin' (l?o) 



1 - £d + [cos2 (^?o) /.J, + sin' (^?o) hy] 



^4 

2no«/^') _ 2 2 ~ + ■ Af [cos' (i?o) /^/j/z + sin' {'do) hyy] 



1 - £d + [cos' (^o) hy + sin' (^?o) /.y] 



2rio^(^) _ 2 ~ ^-d) + 1^^ ■ § [3 cos' (?9o) /y^^ + sin' {do) hyz] 

' 1 - + [cos2 (i?o) Izy + sin' (l?o) Ly] 



_ 2 fgpsin {2do) ^§Jxy, 

^* ' ' 1 - £d + ^ ^ [cos' (^?o) + Sin' {do) Ixy] 



(C2) 



(C3) 
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-^^ 1 - £d + [cos2 (z?o) + sin2 (t?o) 4^] 



2no^(^) ^ 2 ~ + ■ ^ ^f^^ + ^^"^^ ^^2^^] 



2no5^ _ 2 3 (1 - £d) + I^D ■ [5 cos^ (??o) /..z + sin^ (i?o) hzz] 



'XI i-eD + ^§ [cos2 (^o) Izy + siii^ (^?o) 4,] 



2no5^ _ 2 |£Dsin(2^?o)^|/... 



1 - £d + [cos2 (^?o) + sin^ (^?o) 4y] 



2W^) _ 2f,^>^l\ len-sm{2^o)^^h.z 



— / / ^ J « 



2W^)^ 2 A An |£^-sin(2t?o)^|j... 



' 'V A2 ; 1 _ + 3|o ^ [C0S2 (^?o) + Sin^ (^0) hy] 



(C4) 



(C5) 



2no^(^) _ 2 /^A,2 ^ A,2^^ (1 - £d) + f^D ■ ^ [cos^ (t^q) Ixzz + sin^ (t?o) /^x.] 



X A^y 1 - £d + [C0S2 (7?o) /.J/ + sin^ (^0) hy] 
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(C6) 

_ ^2 / (1 - sp) + fgfll [3 cos^ (^o) ly,, + sin^ (^p) 7,,,] 

''''' ~ V 1 - £,5 + [cos2 (^?o) + sin^ (^o) 



2rio5( 



(^) 



-a 



1 + ^ 
A2 



kz5sin(2^o)4-i/ 



1 - + [cos2 (^o) hy + sin^ (^?o) 



(C7) 



\ ^XV.VZ — '-^n, \ J- + 



l-SD + ^t^ [cos2 (^?o) + sin^ (^?o) /.J,] 



Here, the quantities lab and labc denote (scaled) double- and triple index integrals, as ex- 
plained in flAlOp . 



Appendix D: Coupled Monopole-Quadrupole Modes of Density Oscillations 

For completeness, we discuss here the coupled small amplitude monopole-quadrupole 
oscillations of density for a BEC confined in a harmonic trap with cylindrical (uniaxial) 
symmetry, restricting to the case of zero dipole-dipole interaction, Sd = Setting ujz 7^ 
ujy = Ux = in the eigenvalue problem (11161) we easily find analytical expressions for three 
eigenmodes. First 

^^11,2 = V2uj^ (Dl) 
1 

-1 




PXX \P'^j_y2 

Pyy (^^i^-yi 
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It follows directly from (IHTj) that this eigenmode corresponds for all anisotropy ratios to a 
density fluctuation 6nQ (r, t) with pure (ia;2_j^2-symmetry: 



9 9 

5nn (r, t) = 2no cos (fit + 6n) ^ 



A 



(0) 



The second and third eigenmodes 'paa (^^±^ j with eigenfrequencies flf.'^ form a doublet 
consisting of a combination of basis elements with s-wave and dz^-wave symmetry. We 
obtain as a function of the anisotropy ratio = ^ the following exact results for the 
eigenfrequencies and the eigenvectors: 



+ 

(0) 



4-3i/^+Vl6--16v^+9i? 



4~3i^^+Vl6-16i/^+9i^ 



(D2) 



For z/ — > oo this mode becomes quasi one- dimensional 



V » I 



Pyy 

Pzz (^fi 



(.))■ 




1 








?') 




1 








,0. 







+ ... 
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, while for — > it is quasi two-dimensional: 



= 2a;^ ( 1 + 1^.^ + .... 



Pxx 


(fi™ ) ■ 




2 

iy2 


- ^ + 


Pyy 


H"' ) 




2 


4 ••• 


Pzz 


hf 






1 



For an anisotropy ratio ~ 1 (slightly deformed sphere) the associated density fluctuation 
Suq (r, t) is of the breather type, i.e. a strongly weighted isotropic s-wave part is combined 
with only a small admixture of quadrupolar dz^-wsive symmetry: 



(0) 



5nQ (r, t) = 2no cos (Qt + 5n) 



/ 4-3i/^+A/16-16i^^+9i/^' j_ l\ 

y ^2) 



' X ' ' y 

WT 



I f 4-3i/^+\/16-16i^^+9i/'^" _|_ r1 



( 



4i/2 



+ 







(D3) 



The other eigenmode of the doublet is characterized by: 



^(0) /' 4 + 3i/2 - V16 - 16z^2 + 9^,4 



(D4) 



P.. ) _ 



4i^2 



4-3i/-- V'l(-i-lK>y-'+9i? 
4iy2 
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For z/ — )■ oo this mode behaves asymptotically like 



10 / 1 



For z/ — )■ we find 



Pxx 




Pyy 




Pzz 





n 



(0) 



2 ^ 3i/2 ^ •• 



2 + 31/2 + •• 



Pxx 




Pyy 




Pzz 





32 ^ 



32*^ ^ 



1 

For an anisotropy ratio ^ 1 (slightly deformed sphere) this mode describes a density 
fiuctuation with a strongly weighted d^^- wave part and only a small admixture of isotropic 
s-wave symmetry: 



(D5) 



6nn (r, t) = 2no cos (fit + 5^; 



2i/2 2 y i^my 



+ 



4-31/2- Vl6-16i/2+9i/'^ 



4i/^ 



+ 



' z 

A? 



4-31/2- Vl6-16i/2+9i/"4 I 1 
4i/2 2 

The derived frequencies for the coupled monopole-quadrupole oscillations of a BEC without 
dipole-dipole interaction, i.e. e^i = 0, that is confined inside a harmonic trap with uniaxial 
(cylindrical) symmetry, coincide with well known results first derived by Stringari ll|] using 
a different method, that enabled him also to derive all the higher lying frequencies. 
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